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A dynamic crack tip equation of motion is proposed based on the autonomy of the near-tip 
nonlinear zone of scale £„i, symmetry principles, causality and scaling arguments. Causality implies 
that the asymptotic linear-elastic fields at time t are determined by the crack path at a retarded 
time t— Td, where the delay time scales with the ratio of £ n i and the typical wave speed c n i within 
the nonlinear zone. The resulting equation is shown to agree with known results in the quasi-static 
regime. As a first application in the fully dynamic regime, an approximate analysis predicts a 
high-speed oscillatory instability whose characteristic scale is determined by £ n i . This prediction is 
corroborated by experimental results, demonstrating the emergence of crack tip inertia-like effects. 



Introduction. — Fundamental puzzles in the dynamic 
fracture of brittle materials remain unresolved mainly 
due to the lack of a well-established equation of motion 
for a crack's tip. The study of dynamic fracture has fo- 
cused on the central idea of linear-elastic energy flowing 
into the crack tip nonlinear and dissipative zone [D, 0] • 
While this approach is successful in determining the crack 
growth rate when its path is known a priori, it is fun- 
damentally deficient in the general and most interesting 
case in which the crack's path is selected dynamically 
(e.g. instabilities [3, 13]), without being supplemented 
with a path selection rule. This long-standing problem 
hampers the development of a predictive and complete 
theory of dynamic fracture. 

In this Letter we propose a dynamic crack tip equa- 
tion of motion for isotropic materials under plane defor- 
mation, based on rather general physical considerations. 
A basic starting point is the concept of the autonomy 
of the crack tip nonlinear zone in the canonical theory 
of fracture, linear-elastic fracture mechanics (LEFM) [l| . 
The idea is that the mechanical state within the small 
near-tip nonlinear zone of scale £ n i ( "inner problem" ) , 
where LEFM breaks down, is uniquely determined by 
the asymptotic linear-elastic fields surrounding it ( "outer 
problem"), but is otherwise independent of the applied 
loadings and the geometric configuration (e.g. crack 
path) in a given problem. Therefore, the near-tip non- 
linear zone is coupled to the applied loadings and the 
geometric configuration through the asymptotic linear- 
elastic fields and the crack tip itself evolves according to 
the dynamics within the near-tip nonlinear zone. What 
role then plays the near-tip nonlinear zone in determining 
the path selected by a crack tip? 

The main idea of this Letter is that in the presence of a 
finite nonlinear near-tip zone, causality implies that the 
asymptotic linear-elastic fields at a given time t, which 
control the crack tip motion at that time, are determined 
by the crack path at a retarded time t—Td, where the 
delay time scales with the ratio of £ n i and the typical 
wave speed c n i within the nonlinear zone. That is, we 
propose that the only essential properties of the near-tip 
nonlinear zone are its size £ n i and its typical wave speed 



c„i, and that these appear in a macroscopic continuum 
theory mainly through their causal effect. This physical 
effect is missing in LEFM since its basic tenet is that 
the nonlinear near-tip zone acts as an energy sink, but 
otherwise £ n i — > can be assumed. 

The mathematical formulation of these ideas leads to 
a simple, continuum level, dynamic equation of motion 
for a crack's tip. This equation agrees with the well- 
known "principle of local symmetry" Q in the quasi- 
static limit, but is shown to have novel implications in 
the fully dynamic regime. As an example, we perform 
an approximate linear stability analysis of rapid mode I 
cracks. It predicts the existence of a spontaneous sym- 
metry breaking, high-speed oscillatory instability whose 
characteristic scale is determined by £ n i. Using the re- 
cently developed weakly nonlinear dynamic fracture the- 
ory [H, @, 0] to estimate £ n i , this prediction is shown to 
agree well with recent experiments [3[ ■ These results ex- 
plicitly demonstrate the importance of the lengthscale £„i 
at high propagation speeds, as well as the emergence of 
crack tip inertia-like effects. 

Crack tip dynamics. — The mathematical formulation 
of the ideas described above follows in three steps. Con- 
sider a crack under plane deformation conditions, whose 
path is described by r Up (t) and whose tip is surrounded 
by a small nonlinear zone £ n i, see Fig. [TJ Note that £ n i 
is a dynamic quantity that depends, for example, on the 
crack speed, v = \r Up \ (cf. Figs. 1 in [5] and [6]). Out- 
side the nonlinear zone, in an annulus whose width is a 
few times t n iiv), the stress tensor <x is properly described 
by the linear-elastic universal, asymptotic fields [H 
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dinates systems, and i,j run over the polar coordinates 
r and ip. Kjji are the mode I and II stress intensity 
factors respectively, and H I,n (ip,v) are known tensorial 
functions [l| . The concept of autonomy implies that the 
near-tip nonlinear zone is coupled to the large scales only 
through the stress intensity factors Ki t u, which drive 
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The third step in the derivation follows from the obser- 
vation that 0(t) is defined at the crack tip, while Knit) 
is defined a distance l n i away from it. Therefore, causal- 
ity implies that Kn at time t cannot be affected by the 
crack faces created in the time interval [t— T<j,t], with 



(6) 



FIG. 1: A crack with a small nonlinear zone of scale i n i. The 
angle 6 the crack makes with respect to the x-direction is 
shown at two times separated by a delay Ti- 



the crack's tip growth. Therefore, our aim is to derive an 
equation of motion based on these quantities. 

As a first step in this derivation, we follow closely the 
reasoning of [8j]. Denote by t and n the tangent and 
normal unit vectors at the crack tip, respectively (see 
Fig. [T|) and consider the discrete symmetry operation R n 
that transforms n— * — n. Under this symmetry operation 
the relevant quantities of the asymptotic LEFM fields of 
Eq. transform as follows: (i) Kj — > Kj (ii) Kn — > 
—Kn (iii) v — > v. Assuming material isotropy, one can 
write down the most general first order equations that are 
invariant under R n . The first equation is just a kinematic 
relation for the rate of crack tip growth 



d t r tip = v(K I (t),K II (t))t . 
The second one describes crack tip rotation 
d t i oc K H (t)h , 
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where the proportionality coefficient is a true scalar. 

The second step in the derivation amounts to estimat- 
ing the proportionality coefficient in Eq. ([3]) by dimen- 
sional considerations. The existence of a lengthscale l n i, 
the crack speed v and a typical propagation stress inten- 
sity factor K c , imply that 



dtt 



v K H (t) 

'inl K c 



n 



(4) 



K c can be related to the fracture energy T(v), v(Ki, Kn) 
is determined by standard energy balance considerations, 
i.e. the generalized Griffith criterion and the minus 
sign ensures that cracks rotate in the proper direction 
in the presence of mode II fields. Our scaling approach 
assumes that all other material-specific properties of the 
nonlinear zone appear as a pre-factor of order unity in 
Eq. (0| . Equation (0| can be rewritten in terms of the 
angle 9 that the unit tangent t makes with the x-axis as 



Here c n i is the typical wave speed within the nonlinear 
zone, possibly of the order of the linear elastic wave speed 
c s , but not necessarily so. We note that using a single 
delay time T d is certainly a simplification of more compli- 
cated dynamics, but this simplified scaling assumption is 
expected to capture the essence of the physics involved. 

To formulate this idea precisely, we should express the 
physical Kn at time t in terms of the mathematical ICn at 
a retarded time t— Td, taking into account the fact that 
the latter is defined with respect to a coordinates system 
rotated by 8(t — Td), while the former with respect to a 
coordinates system rotated by 9{t), see Fig. [TJ K,n(t—Td) 
is obtained from a pure LEFM problem with a crack path 
corresponding to t—T d . Therefore, we have 



K n (t) ~ )C I (t-T d )Z I r(p (<p = 6(t)-9(t-T d ),v) 

+ lC II (t-Td)X I r l( l p = 6(t)-6(t-Td),v), (7) 



where S^j (ip, v) were defined in Eq. JT]) . The effect of 
the translation of the tip during the time interval Td may 
be non-negligible in general, but for the present purposes 
we neglect it. By combining Eqs. ©-J?]), we obtain the 
proposed dynamic equation of motion for the crack tip. 

The quasi-static limit. — It is clear that in the quasi- 
static limit, v <C c n i , the effect of the delay time Td is negli- 
gible and Eq. predicts Kn=Kn (use 0(t)= 0(t-T d ), 
E T r(p (0,v) = and S^(0,u) = 1 [1]). Moreover, in this 
limit the crack tip has enough time to accommodate the 
presence of Kn and thus (excluding crack initiation un- 
der imposed finite mode-mixity) , we expect a predomi- 
nantly mode I propagation with K c ~ Kj. Scaling out 
v by introducing the arc-length parametrization s of the 
crack path, ds = vdt, Eq. ([5]) becomes 
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The left hand side of Eq. ([5]) has the dimensions of an 
inverse length. In the quasi-static limit, where the effect 
of the delay time Td is negligible, this lengthscale cannot 
be determined by i n i. Therefore, it must be determined 
by a macroscopic lengthscale L that characterizes the 
sample geometry, i.e. \d s 0\ ~ L" 1 . The latter result is 
corroborated by various experimental observations, see 
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for example [ll| , and is consistent with Eq. @ if Kjj / Kj 
is of the order of t n i/L. Therefore, we have 



velocity-dependent functions that appear in it are 



\d 8 6\~L- 



Kn/Kj-Oiint/L)^! 



(9) 

The last result is the celebrated "principle of local sym- 
metry" 0], which is inevitable from the LEFM perspec- 
tive in which l n i does not exist. This principle, coupled 
to the Griffith criterion [l|, provides an excellent quan- 
titative description of quasi-static crack propagation in 
brittle materials 

HE El. In fact > in ^ 

was shown 

that Eqs. ((8]) and the "principle of local symmetry" gen- 
erate indistinguishable predictions. 

High-speed oscillatory instability. — In situations in 
which v is of the order of c n i, the delay time r d may 
be important, possibly giving rise to novel physical ef- 
fects. To explore this possibility we use Eq. ([5|), with 
Eqs. j6|)-([7]), to study the linear stability of rapid mode I 
cracks propagating nearly steadily in a large body of size 
r^L 2 with negligible wave interactions with the bound- 
aries. These conditions can be easily met when an initial 
seed crack accelerates quickly to the center of the body 
[H. Thus, the crack length I is 0(L). 

Consider then a configuration C e that results from 
a small time-dependent perturbation of the crack path 
C e = {(x,y): —L/2 < x < vt, y = eip(x)}. Here ijj{x) is 
a smooth dimensionless function that defines transverse 
perturbations, x is the propagation direction and y is the 
loading direction in which a uniform tensile stress a 00 is 
applied. The dimensional amplitude e is much smaller 
than any other lengthscale in the problem, ensuring that 
both the speed and the path are only slightly perturbed. 
Note that we choose I = L/2 at t = 0. By symmetry, ICi 
and K,n in the perturbed configuration have the forms 
/C/~/c[ 0) +O(e 2 ) and K, II ~K { ^+0{e i ), where the super- 
scripts denote orders in e. The crack speed v is quadratic 
in the stress intensity factors, implying that the mode II 
contribution to v is negligible to 0(e). 

We are interested in situations in which the typical 
scale A of transverse perturbations ip( x ) (soon to be iden- 
tified with a wavelength) satisfies e<A<(, L. The stress 
intensity factors in the mathematical LEFM problem for 
the crack configuration C e were calculated in the gen- 
eral three-dimensional case in 14j and applied to a two- 
dimensional situation in [la ] . The result which is relevant 
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^(-(v) - vQ(v)q n (v))e(t-T d ) 
eB{il>) . (10) 



In this expression, as compared to the one appeared in 
[lit , we omitted a term of order a°°e/ \fl that is negligible 
with respect to 9K { p ~ er°°e\/I/A (e < A < I), where 
ta.n(9) = ed x ip(vt)~9. 

The first term on the right-hand-side (RHS) of Eq. 
([TrJ|) is local, i.e. it depends on the angle 9, and the 
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Q(v) = [2a d (a s -a d ){l+a 2 )] /R(v), 
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and qn(v) is given explicitly in [l5|, |lg]. R(v) is the 
Rayleigh function [3, El, [r| , a 2 d = l—v 2 /c 2 d and c s and 
c d are the shear and dilatational wave speeds, respec- 
tively. B{4>) on the RHS of Eq. (jTUJ) is a non-local linear 
functional (involving a spatio-temporal convolution). It 
contains a contribution that is negligible compared to the 
local term and a contribution that is strictly divergent 

nana. 

The latter was proposed to be regularized in 
terms of generalized functions [la ] and in that case may 
be of the same order of magnitude as the local term [l6| . 
We do not, however, consider the possibly relevant part 
of B(ip) here and neglect it in Eq. (fT0| . We note that 
in the quasi-static limit B(xj)) has a destabilizing role in 
a crack path stability analysis [3, [13] , such that its omis- 
sion here may in fact promote stability. Below we show 
that the local terms are sufficient to induce an instability. 

The resulting equation is substituted in the RHS of Eq. 
([7|) , which is then expanded to first order in 9 (with K, d = 
/Cj°^ independent of t) and the result substituted in the 
RHS of Eq. {5J . This yields the following approximation 
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[E(v)6{t) - vG(v)qn(v)6(t - r d )\ , 



'.K^p was used. 
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where K c 

In order to study the linear stability of Eq. (Q3 
consider linear modes 9(t) — ae luJt , with a ~ e/A and 
5R(w) — 2ttv/\. Here A explicitly denotes the wavelength 
of spatial perturbations. Substituting the linear modes 
into Eq. (p~2|) . we obtain 

(13) 



icniuj = -(3v [S(u) - vQ(v)q u {v)e~ 



where uj = toT d is a dimensionless complex frequency and 
(3 = T d Cni / f-ni is 0(1). The solution Q(v/c s ), which de- 
pends rather weakly on (3, c n i/c s and c d /c s , is shown in 
Fig. [21 The major result is that 3(w) becomes negative 
at a critical speed v c ~ 0.77c s (for /3 = 0.5, c d = 2c s and 
Cni = Cs- The latter is explicitly derived in [lij for the 
material used in 0) with sft(w) t^O. This implies an oscil- 
latory instability at v c . At v c we have 5ft (a) c ) ~ 1.5, which 
implies 



2<irv c pe n i(v c ) 
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The last estimate is obtained using the numbers used in 
Fig. [2l To conclude, the dynamic crack tip equation of 
motion predicts an oscillatory instability at a high critical 
speed v c /c s (weakly material dependent), with an initial 
wavelength A that scales with the linear size of the near- 
tip nonlinear zone £ n i{v c )- 
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FIG. 2: Q(u>) (left) and SR(w) (right) as a function of v/c s 
for f3 — 0.5, c n i = c s and = 2c s 16]. A linear oscillatory 
instability is predicted at t> c ~0.77c s , for which Sr(w) becomes 
negative with 5R(w) 7^ 0. A multivalued function corresponds 
to multiple solutions of Eq. (|13|l . The symmetric 5R(cj) < 
branch is not shown. 



Comparison to experiments. — A high-speed oscilla- 
tory instability in a brittle material was observed ex- 
perimentally in 0]. The reported critical speed was 
i> c ~0.87c\j, which reasonably agrees with our theoretical 
prediction of v c ~ 0.77c s , especially in light of the vari- 
ous scaling and simplifying assumptions adopted above. 
In order to compare the theoretical prediction for the 
wavelength to the measured one, we need an estimate of 
^ni(vc) in Eq. (fT4"|) . The recently developed weakly non- 
linear dynamic fracture theory 0, @] predicts nonlinear 
corrections to LEFM, cf. Eq. (fTJ), on a dynamic scale 
£ni(v) — h(v)G(v) I ji. h(v) is a material-dependent func- 
tion, G(v) is the LEFM energy release rate that balances 
the fracture energy T(v) [l| and fi is the shear modulus. 
Using the numbers reported in Fig. 1 of [6[ , we estimate 
4i(w c )~4mm [16l | . which upon substitution in Eq. (I14p 
yields A ~ 6.5mm. 

The saturated (nonlinear) wavelength X final, which 
cannot be computed within our linearized analysis, is re- 
ported in Fig. 4b of [3| and shows a systematic variation 
with a°° . The initial wavelength A^ (i.e. at the onset 
of oscillations), however, shows no systematic variation 
with er°°, attaining the value Aj n , = 7.5±2.3mm [17j . This 
experimental result is in agreement with our predicted A, 
which pertains to the linearized dynamics near the onset 
of instability. Furthermore, Ref. Q reports that when 
the shear wave speed of the material used in the exper- 
iments was increased by a factor of 3, no appreciable 
change in the reduced critical speed v c /c s was observed, 
in excellent agreement with the present predictions. Fi- 
nally, we note that e~10 _1 mni (e corresponds to Ai n i in 
the notation of I3j), A ~ 7.5mm and L~ 10 2 mm in the 
experiments of [3j, are all consistent with the required 
scales separation £<A<L. To conclude, the theoretical 
predictions of the proposed theory are in agreement with 
the available experimental data. 

Summary. — We have proposed a dynamic equation 



of motion for crack tips propagating in isotropic mate- 
rials under plane deformation. A crucial ingredient in 
the derivation is the idea that a finite near-tip nonlinear 
zone £ n i introduces a delay time Td between the driving 
stress intensity factors and the crack tip itself. The re- 
sulting equation agrees with known quasi-static results 
and predicts a high-speed oscillatory instability with a 
characteristic scale £ n i, in the fully dynamic regime. The 
latter result shows explicitly that the non-geometric, dy- 
namic lengthscale £ n i(v) directly influences the crack's 
dynamics at high speeds, cf. [xo| . Furthermore, the exis- 
tence of a finite i n \ (y) induces crack tip inertia-like effects 
associated with the delay time Td- All of these effects are 
missing in LEFM 0,0]. 

When the present predictions are supplemented with 
the predictions of the recently developed weakly nonlin- 
ear fracture theory for £ n i(y), a satisfactory agreement 
with experimental results for the onset of a high-speed os- 
cillatory instability is obtained. Therefore, extending the 
present ideas to three-dimensional situations may open 
the way to understanding the side-branching instability, 
which also involves a poorly understood lengthscale Q ■ 
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